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ABSTRACT 

We present the evolution of the rest-frame optical luminosity density, j'^)^^^, of the integrated rest- 
frame optical color, and of the stellar mass density, p*, for a sample of X^-band selected galaxies 
in the HDF-S. We derived in the rest-frame U, B, and V^-bands and found that j'^^^^ increases 
by a factor of 1.9 ± 0.4, 2.9 ± 0.6, and 4.9 ± 1.0 in the V, B, and U rest-frame bands respectively 
between a redshift of 0.1 and 3.2. We derived the luminosity weighted mean cosmic {U — B)rest and 
(B — V)rest colors as a function of redshift. The colors bluen almost monotonically with increasing 
redshift; at z = 0.1, the {U — B)rest and {B — V)rest colors are 0.16 and 0.75 respectively, while at 
z = 2.8 they are -0.39 and 0.29 respectively. We derived the luminosity weighted mean M/Ly using 
the correlation between ([/ — V)rest and logioM/L^ which exists for a range in smooth SFHs and 
moderate extinctions. We have shown that the mean of individual M/Ly estimates can overpredict 
the true value by ~ 70% while our method overpredicts the true values by only 35%. We find that 
the universe at z 3 had 10 times lower stellar mass density than it does today in galaxies with 
1V> 1.4 X 10^" h^Q^ Lq. 50% of the stellar mass of the universe was formed by z ~ 1 — 1.5. The rate 
of increase in with decreasing redshift is similar to but above that for independent estimates from 
the HDF-N, but is slightly less than that predicted by the integral of the SFR(z) curve. 
Subject headings: Evolution — galaxies: formation — galaxies: high redshift — galaxies: stellar 
content — galaxies: galaxies 



1. INTRODUCTION 

A primary goal of galaxy evolution studies is to eluci- 
date how the stellar content of the present universe was 
assembled over time. Enormous progress has been made 
in this field over the past decade, driven by advances over 
three different redshift ranges. Large scale redshift sur- 
veys with median redshifts of z ~ 0.1 such as the Sloan 
Digital Sky Survey (SDSS; York et al. 2000) and the 2dF 
Galaxy Redshift Survey (2dFGRS; CoUess et al. 2001), 
coupled with the near infrared (NIR) photometry from 
the 2 Micron AU Sky Survey (2MASS; Skrutskie et al. 
1997), have recently been able to assemble the complete 
samples, with significant co-moving volumes, necessary 
to establish crucial local reference points for the local lu- 
minosity function (e.g. Folkes et al. 1999; Blanton et al. 
2001; Norberg et al. 2002; Blanton et al. 2003c) and the 
local stellar mass function of galaxies l)Bell et alJl2003at 
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ICole et al.ll2nol . 

At z < 1, the pioneering study of galaxy evolution was 
the Canada France Redshift Survey (CFRS; Lilly et al. 
1996). The strength of this survey lay not only in the 
large numbers of galaxies with confirmed spectroscopic 
redshifts, but also in the /-band selection, which enabled 
galaxies at z < 1 to be selected in the rest-frame optical, 
the same way in which galaxies are selected in the local 
universe. 

At high redshifts the field was revolutionized by the 
identification, and subsequent detailed follow-up, of a 
large populat i on of star-forming galaxies at z > 2 
ijSteidel et alJ Il99(j) . These Lyman Break Galaxies 
(LBGs) are identified by the signature of the redshifted 
break in the far UV continuum caused by intervening 
and intrinsic neutral hydrogen absorption. There are 
over 1000 spectroscopically confirmed LBGs at z > 2, 
together with the analogous U-dropout galaxies identi- 
fied using Hubble Space Telescope (HST) filters. The 
individual properties of LBGs have been studied in great 
detail. Estimates for their star formation rates (SFRs), 
extinctions, ages, and stellar masses have been estimated 
by modeling the broad band fiuxes (Sawicki & Yee 1998; 
hereafter SY98; Papovich et al. 2001; hereafter POl; 
Shapleyetal. 2001). Independent measures of their kine- 
matic masses, metallicities, SFRs, and initial mass func- 
tions (IMFs) have been de termined using rest-frame UV 
and optical spcctro scopv l|Er b et al.' 120031 : iPettini et all 
12000.. .2001. 2002: Shaolev et al.,..200ljl . 

Despite these advances, it has proven difficult to rec- 
oncile the ages, SFRs, and stellar masses of individual 
galaxies at different redshifts within a single galaxy for- 
mation scenario. Low redshift studies of the fundamental 
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plane indicate that the stars in elliptical galaxies must 
have been formed by z > 2 (e.g., van Dokkum et al. 
2001) and observations of evolved galaxies at 1 < z < 2 
indicate that the present population of elliptical galaxies 
was already in place &t z > 2.5 (e.g., Benitez et al. 1999; 
Cimatti et al. 2002; but see Zepf 1997). In contrast, 
studies of star forming Lyman Break Galaxies spectro- 
scopically confirmed to lie at z > 2 (LBGs; Steidel et al. 
1996, 1999) claim that LBGs are uniformly very young 
and a factor of 10 less massive than present day L^,, galax- 
ies (e.g., SY98; POl; Shapley et al. 2001). 

An alternative method of tracking the build-up of the 
cosmic stellar mass is to measure the total emissivity 
of all relatively unobscurcd stars in the universe, thus 
effectively making a luminosity weighted mean of the 
galaxy population. This can be partly accomplished by 
measuring the evolution in the global luminosity density 
j{z) from galaxy redshift surveys. Early studies at in- 
termediate redshift have shown that the rest-frame UV 
and B-band j{z) are steeply increasing out to z ~ 1 
(e.g., Lilly et al. 1996; Fried et al. 2001). iWolf et all 
(|200^ has recently measured j{z) at < z < 1.2 from 
the COMBO-17 survey using ~25,000 galaxies with red- 
shifts accurate to ~ 0.03 and a total area of 0.78 de- 
grees. A t rest-frame 2800A these measurements confirm 
those of iLillv et alJ l)1996j) but do not support claims 
for a shallower incre ase with redshift which goes like 
(1-l- z)^-^ as claime d bvlCowie. Songaila. fc Barged l|T999D 
and iWilson et aO ll2002f). On the oth er hand, the B- 
band evolution from iWolf et alJ 1)200 3'^ is only a factor 
of ~ 1.6 between < z < 1, considerabl y shallower 
than the factor of ~ 3.75 increase seen by iLillv et al] 
()1996j) . At z > 2 measurements of the rest- frame UV 
j{z) have been made using the optically selected LBG 
samples (e.g., Madau et al. 1996; Sawicki, Lin, & Yee 
1997; Steidel et al. 1999; Poh et al. 2001) and NIR se- 
lected samples iKashikawa et all 120031 iPoH et alJl2fl(i^ 
iThompsonl l200,'ll and, with modest extinction correc- 
tions, the most recent estimates generically yield rest- 
frame UV j{z) curves which, at z > 2, are approximately 
flat out to z '--^ 6 (cf. Lanzetta et al. 2002). Dickinson et 
al. (2003; hereafter D03) have used deep NIR data from 
NICMOS in the HDF-N to measure the rest-frame B- 
band luminosity density out to z ^ 3, finding that it re- 
mained constant to within a factor of ~ 3. By combining 
j{z) measure ments at different rest-frame wavelen gths 
and redshifts. iMadau. Pozzetti. Sz DickinsonI l)1998f) and 
iPei. Fall, fc HauseT ' ifl99'9f) modeled the emission in all 
bands using an assumed global SFH and used it to 
constrain the mean extinction, metallicity, and IMF. 
iBolzonella. Pcllo. & Maccagni ( 2002) measured NIR lu- 
minosity functions in the HDF-N and HDF-S and find 
little evolution in the bright end of the galaxy popu- 
lation and no decline in the rest-fr ame NIR lum i nosity 
density out to z ^ 2. In add ition, iBaldrv et alJ l)2002f l 
and Glazebroo k et al.) l)2003al) have used the mean opti- 
cal cosmic spectrum at z ^ from the 2dFGRS and the 
SDSS respectively to constrain the cosmic star formation 
history. 

Despite the wealth of information obtained from stud- 
ies of the integrated galaxy population, there are major 
difficulties in using these many disparate measurements 
to re-construct the evolution in the stellar mass density. 



First, and perhaps most important, the selection criteria 
for the low and high redshift surveys are usually vastly 
different. At z < 1 galaxies are selected by their rest- 
frame optical light. At z > 2, however, the dearth of 
deep, wide- field NIR imaging has forced galaxy selec- 
tion by the rest-frame UV light. Observations in the 
rest-frame UV are much more sensitive to the presence 
of young stars and extinction than observations in the 
rest-frame optical. Second, state-of-the-art deep surveys 
have only been performed in small fields and the effects 
of field-to-field variance at faint magnitudes, and in the 
rest-frame optical, are not well understood. 

In the face of field-to-field variance, the globally av- 
eraged rest-frame color may be a more robust charac- 
terization of the galaxy population than either the lu- 
minosity density or the mass density because it is, to 
the first order, insensitive to the exact density normal- 
ization. At the same time, it encodes information about 
the dust obscuration, metallicity, and SFH of the cosmic 
stellar population. It therefore provides an important 
constraint on galaxy formation models which may be re- 
liably determined from relatively small fields. 

To track consistently the globally averaged evolution 
of the galaxies which dominate the stellar mass budget 
of the universe - as opposed to the UV luminosity bud- 
get - over a large redshift range a different strategy than 
UV selection must be adopted. It is not only desirable 
to measure j(z) in a constant rest-frame optical band- 
pass, but it is also necessary that galaxies be selected by 
light redward of the Balmer/4000A break, where the light 
from older stars contributes significantly to the SED. 
To accomplish this, we obtained ultra- deep NIR imag- 
ing o f the WFPC2 field of the HDF-S l|Ca,Rerta,no et al J 
120001) with the Infrared Spectrograph And Array Camera 
(ISAAC; Moorwood et al. 1997) at the Very Large Tele- 
scope (VLT) as part of the Faint Infrared Extragalactic 
Survey (FIRES; Franx et al. 2000). The FIRES data 
on the HDF-S, detailed in Labbe et al. (2003; hereafter 
L03), provide us with the deepest ground-based Js and 
H data and the overall deepest i^^-band data in any 
field allowing us to reach rest-frame optical luminosities 
in the V-band of ~ 0.6 at z ~ 3. First results 

using a smaller set of the data were presented in Rud- 
nick et al. (2001; hereafter ROl). The second FIRES 
field, centered on the z = 0.83 cluster MS1054-03, has 
^ 1 magnitude less depth bu t ~ 5 times greater area 
l|Forster Schreiber et a,].ll200l. 

In the present work we will draw on photometric red- 
shift estimates, Zphot ior the X^-band selected sample in 
the HDF-S (ROl; L03), and on the observed SEDs, to 
derived rest-frame optical luminosities L^°''* for a sam- 
ple of galaxies selected by light redder than the rest- 
frame optical out to z ~ 3. In § |21 we describe the 
observations, data reduction, and the construction of 
a i^s-band selected catalog with 0.3 — 2.2/im photom- 
etry, which selects galaxies at z < 4 by light redward 
of the 4000A break. In § we describe our photomet- 
ric redshift technique, how we estimate the associated 
uncertainties in Zphot, and how we measure L^'^''* for 
our galaxies. In § ^ we use our measures of 17^^^ for 
the individual galaxies to derive the mean cosmic lumi- 
nosity density, f^^*" and the cosmic color and then use 
these to measure the stellar mass density p* as a func- 
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tion of cosmic time. We discuss our results in § |5l and 
summarize in § Throughout this paper we assume 
r^M = 0.3, riA = 0.7, and Ho = 70 hyo km s-^Mpc"! 
unless explicitly stated otherwise. 

2. DATA 

A complete description of the FIRES observations, re- 
duction procedures, and the construction of photometric 
catalogs is presented in detail in LOS; we outline the im- 
portant steps below. 

Objects were detected in the j^Ts-ba nd image with ver- 
sion 2.2.2 of the SExtractor software l|Bertin fc Arnout^ 
Il996j) . For consistent photometry between the space and 
ground-based data, all images were then convolved to 
0'.'48, the seeing in our worst NIR band. Photometry 
was then performed in the [/300, S450, Vgoe, hu, Js, H, 
and Ks-hand images using specially tailored isophotal 
apertures defined from the detection image. In addition, 
a measurement of the total flux in the A's-band, -f^^*°AB' 
was obtained using an aperture based on the SExtractor 
AUTO aperture^. Our effective area is 4.74 square ar- 
cminutes, including only areas of the chip which were well 
exposed. All magnitudes are quoted in the Vega system 
unless specifically noted otherwise. Our adopted conver- 
sions from Vega system to the AB system are Js,vega = 
AR - 0.90, H„p„n. = Hah - 1.38, and Kg^vega = Ks^ab - 
1.86 lIBessell fc BrettllT988l) . 

3. MEASURING PHOTOMETRIC REDSHIFTS AND 
REST-FRAME LUMINOSITIES 

3.1. Photometric Redshift Technique 

We estimated Zphot from the broad-band SED using 
the method described in ROl, which attempts to fit the 
observed SED with a linear combination of redshifted 
galaxy templates. We made two modifications to the ROl 
method. First, we added an additional template con- 
structed from a 10 Myr old , single age, solar metallicity 
population with a ISalpeteil l|195,'Tf) initial mass function 
(IMF) based on empirical stellar s pectr a from the 1999 
version of the iBruzual fc Chariot (|199?) stellar popula- 
tion synthesis code. Second, a 5% minimum flux error 
was adopted for all bands to account for the night-to- 
night uncertainty in the derived zeropoints and for tem- 
plate mismatch effects, although in reality both of these 
errors are non-gaussian. 

Using 39 galaxies with reliable FI RES photome- 
try a n d spectroscopy availabl e from 'Cristiani ct alj 
^2^. migonoulon et a,].l Eml. l?M.rookJ.2003M9. 
iVanzella et an()2002|) . and lRudnick et alJ l)2003|) we mea- 
sured the redshift accuracy of our technique to be 
( \zspec - Zphotl / (1 + Zspec) ) = 0.09 for z < 3 . There 
is one galaxy at Zgpec =2.025 with Zphot ~ 0.12 but with 
a very large internal Zp^ot uncertainty. When this ob- 
ject is removed, ( jz^pec - Zphot\ / (1 + ^^spec) ) = 0.05 at 

Zspe.c ^ 1.3. 

For a given galaxy, the photometric redshift probabil- 
ity distribution can be highly non-Gaussian and contain 
multiple minima at vastly different redshifts. An ac- 
curate estimate of the error in Zp^ot must therefore not 

® The reduced images, photometric catalogs, photo- 
metric redshift estimates, and rest-frame luminosities 
are available online through the FIRES homepage at 
http : //www. St rw. leidenuniv.nl/~f ires. 

^ available at http://www.aao.gov.au/hdfs/ 



only contain the two-sided confidence interval in the lo- 
cal minimum, but also reflect the presence of alternate 
redshift solutions. The difficulties of measuring the un- 
certainty in Zphot were discussed in ROl and will not be 
repeated in detail here. To improve on ROl, however, 
we have developed a Monte Carlo method which takes 
into account, on a galaxy- by-galaxy basis, flux errors 
and template mismatch. These uncertainty estimates 
are called Szphot- For a full discussion of this method 
see Appendix 1X1 

Galaxies with Kl°^^ > 25 have such high photometric 
errors that the Zphot estimates can be very uncertain. At 
^s°AB ^ 25, however, objects are detected at better than 
the lO-sigma level and have well measured NIR SEDs, 
important for locating redshifted optical breaks. For this 
reason, we limited our catalog to the 329 objects that 
have i^'°AB < -^5, lie on well exposed sections of the 
chip, and are not identified as stars (see ii3.1.1|l . 

3.1.1. Star Identification 

To identify probable stars in our catalog we did not 
use the profiles measured from the WFPC2-imaging 
because it is difficult to determine the size at faint 
levels. At the same time, we verified that the stellar 
template fitting technique identified all bright un- 
saturated stars in the image. Instead, we compared 
the observed SEDs with those from 135 NextGen 
version 5.0 stellar atmosphere m odels described in 
iHauschildt. Allard. fc BaronI l)1999() and available at 
http : //dilbert .physast .uga.edu/~yeti/mdwarfs .html. 
We used models with log(g) of 5.5 and 6, effective tem- 
peratures ranging from 1600 K to 10,000 K, and 
metallicities of solar and 1/lOth solar. We identified an 
object as a stellar candidate if the raw of the stellar 
fit was lower than that of the best-fit galaxy template 
combination. Four of the stellar candidates from this 
technique (objects 155, 230, 296, and 323) are obviously 
extended and were excluded from the list of stellar 
candidates. Two bright stars (objects 39 and 51) were 
not not identified by this technique because they are 
saturated in the HST images and were added to the list 
by hand. We ended up with a list of 29 stars that had 
^Tab < 25 and lie on well exposed sections of the chip. 
These were excluded from all further analysis. 

3.2. Rest-Frame Luminosities 

To measure the L™*'* of a galaxy one must combine its 
redshift with the observed SED to estimate the intrinsic 
SED. In practice, this requires some assumptions about 
the intrinsic SED. 

In ROl we derived rest-frame luminosities from the 
best-fit combination of spectral templates at Zphot, which 
assumes that the intrinsic SED is well modeled by our 
template set. We know that for many galaxies the best-fit 
template matches the position and strength of the spec- 
tral breaks and the general shape of the SED. There are, 
however, galaxies in our sample which show clear resid- 
uals from the best fit template combination. Even for 
the qualitatively good fits, the model and observed fiux 
points can differ by ~ 10%, corresponding to a ~ 15% 
error in the derived rest-frame color. As we will see in 
ti4.2l such color errors can cause errors of up to a factor 
of 1.5 in the F-band stellar mass-to-light ratio, M/Ly. 
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Here we used a method of estimating l/^^*' which does 
not depend directly on template fits to the data but, 
rather, interpolates directly between the observed bands 
using the templates as a guide. We define our rest- 
frame photometric system in Appendix ^ and explain 
our method for estimating L"*^* in Appendix O 

We plot in Figure^the rest-frame luminosities vs. red- 
shift and enclosed volume for the ^^'"ab < 25 galaxies in 
the FIRES sample. The different symbols represent dif- 
ferent Szphot values and since the derived luminosity is 
tightly coupled to the redshift, we do not independently 
plot U^^^ errorbars. The tracks indicate the L™**' for dif- 
ferent SED types normalized to J^*°ab = 25, while the 
intersection of the tracks in each panel indicates the red- 
shift at which the rest-frame filter passes through our 
i^g-band detection filter. There is a wide range in L^°*'' 
at all redshifts and there are galaxies at z > 2 with U-^^^ 
much in excess of the local L, values. Using the full 
FIRES dataset, we are much more sensitive than in ROl; 
objects at z « 3 with Kl°l^ = 25 have IVw 0.6*Li^°^^i, as 

defined from the z=0.1 sample of Blanton et al. (2003c; 
hereafter B03). As seen in ROl there are many galax- 
ies at z > 2, in all bands, with L™"* >L'°'=^'. ROl 
found 10 galaxies at 2 < z < 3.5 with 1B> lO" hf^f Lq 
and inferred a brightening in the luminosity function of 
~ 1 — 1.3 magnitudes. We confirm their result when 
using the same local luminosity function (Blanto n et al.l 
12001(1 . Although this brightening is biased upwards by 
photometric redshift errors, we find a similar brightening 
of approximately ~ 1 magnitude after correction for this 
effect. As also noticed in ROl, we found a deficit of lumi- 
nous galaxies at 1.5 < z < 2 although this deficit is not 
as pronounced at lower values of L^"*^*. The photomet- 
ric redshifts in the HDF-S, however, are not well tested 
in this regime. To help judge the reality of this deficit 
we compared our photometric redshifts on an objcct-by- 
object basis to those of the Rome group ( Fontana ct al. 
|2000() ^° who derived Zphot estimates for galaxies in the 
HDF-S using much shallower NIR data. We find gen- 
erally good agreement in the Zphot estimates, although 
there is a large scatter at 1.5 < z < 2.0. Both sets of 
photometric redshifts show a deficit in the Zphot distribu- 
tion, although the Rome group's gap is less pronounced 
than ours and is at a slightly lower redshift. In addi- 
tion, we examined the photometric redshift distribution 
of the NIR selected galaxies of D03 in the HDF-N, which 
have very deep NIR data. These galaxies also showed a 
gap in the Zphot distribution at z ^ 1.6. Together these 
results indicate that systematic effects in the Zphot de- 
terminations may be significant at 1.5 < z < 2.0. On 
the other hand, we also derived photometric redshifts for 
a preliminary set of data in the MS1054-03 field of the 
FIRES survey, whose filter set is similar, but which has 
a U instead of U300 filter. In this field, no systematic de- 
pletion of 1.5 < z < 2 galaxies was found. It is therefore 
not clear what role systematic effects play in comparison 
to ficld-to-ficld variations in the true redshift distribu- 
tion over this redshift range. Obtaining spectroscopic 
redshifts at 1.5 < z < 2 is the only way to judge the 
accuracy of the Zphot estimates in this regime. 
We have also split the points up according to 

available at http://www.niporzio.astro.it/HIGHZ/HDF.htiiil 



whether or not they sa t isfied the U-dropout criteria of 
IGiavalisco fc DickinsonI l|200H) which were designed to 
pick unobscured star- forming galaxies at z > 2. As ex- 
pected from the high efficiency of the U-dropout tech- 
nique, we find that only 15% of the 57 classified U- 
dropouts have Zphot < 2. As we will discuss in tl4.1l 
we measured the luminosity density for objects with 
1V> 1.4 X 10^° hfo^ Lq. Above this threshold, there 
are 62 galaxies with 2 < z < 3.2, of which 26 are not 
classified as U-dropouts. These non U-dropouts number 
among the most rest-frame optically luminous galaxies 
in our sample. In fact, the most rest-frame optically lu- 
minous object at z < 3.2 (object 611) is a galaxy which 
fails the U-dropout criteria. 10 of these 26 objects, in- 
cluding object 611, also have J — K > 2.3, a c olor thresh- 
old which has been show n by iFranx et al.l l|2"003il and 
Ivan Dokkum et aP l|2003|) to efficiently select galaxies at 
z > 2. These galaxies are not only luminous but also 
have re d rest-frame optica l colors, implying high M/L* 
values. iFranx et al.l l(200 3^ showed that they likely con- 
tribute significantly (~ 43%) to the stellar mass budget 
at high redshifts. 

3.2.1. Emission Lines 

There will be emission line contamination of the rest- 
frame broad-band luminosities when rest-frame optical 
emission lines contribute significantly to the flux in our 
observed filters. POl estimated the effect of emission lines 
in the NICMOS F160W filter and the filter and found 
that redshifted, rest-frame optical emission lines, whose 
equivalent widths are at the maximum end of those ob- 
served for starburst galaxies (rest-frame equivalent width 
- 200 A), can contribute up to 0.2 magnitudes in the 
NIR filters. In addition, m odels of emission lines from 
iCharlot & Lonefhettil l|2001(l show that emission lines will 
tend to drive the {U — B)rest color to the blue more eas- 
ily than the {B — V)rest color for a large range of models. 
Using the UBV photometry and spectra of nearby galax- 
ies from the Nearby Field Galaxy Survey (NFGS; Jansen 
et al. 2000a; Jansen et al. 2000b) we computed the ac- 
tual correction to the {U — B)rest and {B — V)rest 
colors as a function of (B — V)rest- For the bluest galax- 
ies in {B — V)rest, emission lines bluen the (U — B)rest 
colors by ^ 0.05 and the [B — V)rest colors only by 
< 0.01. Without knowing beforehand the strength of 
emission lines in any of our galaxies, we corrected our 
rest-frame colors based on the results from Jansen et al. 
We ignored the very small correction to the {B — V)rest 
colors and corrected the [U — B)rest colors using the 
equation: 

{U-B),orrected = iU-B)-0M58x{B-V)+0.0656 (1) 

which corresponds to a linear fit to the NFGS data. 
These effects might be greater for objects with strong 
AGN contribution to their fluxes. 

4. THE PROPERTIES OF THE MASSIVE GALAXY 
POPULATION 

In this section we discuss the use of the Zphot and L'^'^'^' 
estimates to derive the integrated properties of the pop- 
ulation, namely the luminosity density, the mean cosmic 
rest-frame color, the stellar mass-to-light ratio M/L*, 
and the stellar mass density p^. As will be described 
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Fig. 1. — The distribution of rest-frame V, B, and {/-band lu- 
minosities as a function of enclosed co-moving volume and ^p^ot 
is shown in figures (a), (b), and (c) respectively for galaxies with 

^tot 



Large symbols have (5zp/io(/(l 



< 0.16 and 



small symbols have <52pf,ot/(l + Zpf^^t) > 0.16. Triangle points 
would be classified as U-dropouts according to the selection of 
IHiavalisco & Dickinson 1 20011). As is expected, most of the galax- 
ies selected as U-dropouts have Zpf^^t > 2. Note, however, the large 
numbers of rest-frame optically luminous galaxies at ^ > 2 which 
would not be selected as U-dropouts. The large stars in each panel 
indicate the value of L',°'^='' from lBlanton et al.l p)03c). In the V- 
band we are sensitive to galaxies at 60% of Ljf^^' even at z ~ 3 
and there are galaxies at Zp^ot ^ 2 with L^''''* > 10^^ hyp Lq. 
The tracks represent the values of L^''''* for our seven template 
spectra normalized at each redshift to K^°p^^ = 25. The specific 
tracks correspond to the E (solid), Sbc (dot), Scd (short dash), Im 
(long dash), SBl (dot-short dash), SB2 (dot-long dash), and lOmy 
(dot) templates. The horizontal dotted line in (a) indicates the lu- 
minosity threshold L^'^'^'^'^ above which we measure the rest-frame 
luminosity density j^^^^ and the vertical dotted lines in each panel 
mark the redshift boundaries of the regions for which we measure 
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below, addressing the integrated properties of the popu- 
lation reduces many of the uncertainties associated with 
modeling individual galaxies and, in the case of the cos- 
mic color, is less sensitive to field-to-field variations. 

4.1. The Luminosity Density 
Using our L'^"*'' estimates from the KI°^q < 25 galaxies 
(see ^3.2|l . we traced the redshift evolution of the rest- 
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frame optically most luminous, and therefore presumably 
most massive, galaxies by measuring the rest-frame lu- 
minosity density of the visible stars associated with 
them. The results are presented in Table ^ ^^'^ plot- 
ted against redshift and elapsed cosmic time in Figure |21 
As our best alternative to a selection by galaxy mass, 
we selected our galaxies in our reddest rest-frame band 
available at z 3, i.e. the y-band. In choosing the z 
and U^^*" regime over which we measured j"''* we wanted 
to push to as high of a redshift as possible with the dou- 
ble constraint that the redshifted rest-frame filter still 
overlapped with the Ks filter and that we were equally 
complete at all considered redshifts. By choosing an IV 



threshold, L*J"°"h 



1.4 X 10^° h 



70 



^0 7 



and a maximum 



redshift of z — 3.2, we could select galaxies down to 
0.6 L'°y ' with constant efficiency regardless of SED type. 
We then divided the range out to z = 3.2 into three bins 
of equal co-moving volume which correspond to the red- 
shift intervals 0-1.6, 1.6-2.41, and 2.41-3.2. 

In a given redshift interval, we estimated jj^'^''* directly 
from the data in two steps. We first added up all the lu- 
minosities of galaxies which satisfied our Ly criteria 
defined above and which had dzphot /(I + Zphot) < 0-16, 
roughly twice the mean disagreement between Zphot and 
Zspec (see il3.1|l . Galaxies rejected by our Szphot cut but 
with lV>Ly"''^, however, contribute to the total lumi- 
nosity although they are not included in this first es- 
timate. Under the assumption that these galaxies are 
drawn from the same luminosity function as those which 
passed the Szphot cut, we computed the total luminosity, 
including the light from the nacc accepted galaxies and 
the light lost from the Urej rejected galaxies as 



X (1 + ::^). 



(2) 



As a test of the underlying assumption of this correction 
we performed a K-S test on the distributions of Ks mag- 
nitudes for the rejected and accepted galaxies in each of 
our three volume bins. In all three redshift bins, the re- 
jected galaxies have KI°^q distributions which are con- 
sistent at the > 90% level with being drawn from the 
same magnitude distribution as the accepted galaxies. 
The total correction per volume bin ranges from 5 — 10% 
in every bin. Our results do not change if wc omit those 
galaxies whose photometric redshift probability distribu- 
tion indicates a secondary minima in chisquared. 
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Uncertainties in the luminosity density were computed 
by bootstrapping from the /^*°ab < 25 subsample. This 
method does not take cosmic variance into account and 
the errors may therefore underestimate the true error, 
which includes field-to-field variance. 

Redshift errors might effect the luminosity density in 
a systematic way, as they produce a large error in the 
measured luminosity. This, combined with a steep lu- 
minosity function can bias the observed luminosities up- 
wards, especially at the bright end. This effect can be 
corrected for in a full determination of the luminosity 
function (e.g., Chen et al. 2003), but for our application 
we estimated the strength of this effect by Monte-Carlo 
simulations. When we used the formal redshift errors we 
obtained a very small bias (6%), if we increase the photo- 
metric redshift errors in the simulation to be minimally 
as large as 0.08 * (1 + z), we still obtain a bias in the 
luminosity density on the order of 10% or less. 

Because we exclude galaxies with faint rest-frame lumi- 
nosities or low apparent magnitudes, and do not correct 
for this incompleteness, our estimates should be regarded 
as lower limits on the total luminosity density. One pos- 
sibility for estimating the total luminosity density would 
be to fit a luminosity function as a function of redshift 
and then integrate it over the whole luminosity range. 
We don't go faint enough at high redshift, however, to 
tightly constrain the faint-end slope a. Because extrap- 
olation of jj^^^ to arbitrarily low luminosities is very de- 
pendent on the value of a, we choose to use this simple 
and direct method instead. Including all galaxies with 
^*°AB < 25 raises the f^'^^ values in the z — — 1.6 red- 
shift bin by 86%, 74%, and 66% in the U, B, and V"-bands 
respectively. Likewise, the j™''* values would increase by 
38%, 35%, and 44% for the z = 1.6 - 2.41 bin and would 
increase by 5%, 5%, and 2% for the z = 2.41 - 3.2 bin, 
again in the B, and 1^-bands respectively. 

The dip in the luminosity density in the second lowest 
redshift bin of the (a) and (b) panels of Figure 12 can be 
traced to the lack of intrinsically luminous galaxies at 
z ~ 1.5 - 2 (iJO ROl). The dip is not noticeable in the 
J7-band because the galaxies at z 2 are brighter with 
respect to the z < 1.6 galaxies in the [/-band than in the 
V or B-band, i.e. they have bluer [U — B)rest colors and 
{U — V)rest colors than galaxies at z < 1.6. This lack 
of rest-frame optically bright galaxies at z '--^ 1.5 — 2 may 
result from systematics in the Zphot estimates, which are 
poorly tested in this regime and where the Lyman break 
has not yet entered the C/300 filter, or may refiect a true 
deficit in the redshift distribution of Ks-hand luminous 
galaxies (see t)3.2|l . 

At z < 1 our survey is limited by its small volume. For 
this reason, we supplement our data with other estimates 
of j^''''* at z < 1. 

We compared our results with those of the SDSS 
as follows. First, we selected SDSS Main sample 
galaxies (Strauss et al. 2002) with redshifts in the 
SDSS Early Data Release (Stoughton et al. 2002) 
in the EDR sample provided by and described by 
iBlanton et al.l ll 2003afl. Us ing the product kcorrect 
vl_16 jBlanton et al.ll2003b|) . for each galaxy we fit an 
optical SED to the o.i.yO. 1^0. 1^0. 1^0.1^ magnitudes, after 
correcting the magnitudes to z = 0. 1 for evolution using 
the results of IBlanton et alJ lj2003(Jl . We projected this 



SED onto the UBV fihers as described by Besseh (1990) 
to obtain absolute magnitudes in the U B V Vega-relative 
system. Using the method described in IBlanton et alJ 
( 2003a) we calculated the maximum volume Knax within 
the EDR over which each galaxy could have been ob- 
served, accounting for the survey completeness map and 
the flux limit as a function of position. l/14nax then rep- 
resents the number density contribution of each galaxy. 
From these results we constructed the number density 
distribution of galaxies as a function of color and abso- 
lute magnitude and the contribution to the uncertain- 
ties in those densities from Poisson statistics. While the 
Poisson errors in the SDSS are negligible, cosmic vari- 
ance does contribute to the uncertainties. For a more 
realistic error estimate, we use the fractional errors on 
the luminosity density from .Blanton et al . (2003c). For 
the SDSS luminosity function, our L^'^^^^ encompasses 
54% of the total light. 

In Figure 121 we also sh ow the .7'^'^''* measu rements from 
the COMBO- 17 survey (Wolf et al."2003). We used a 
catalog with updated redshifts and 29471 galaxies at 
z < 0.9, of which 7441 had W>Lf"'''^ (the J2003 cata- 
log; Wolf, C. private communication). Using this cata- 
log we calculated in an identical way to how it was 
calculated for the FIRES data. We divided the data 
into redshift bins of Az = 0.2 and counted the light 
from all galaxies contained within each bin which had 
lV>Ly The formal 68% confidence limits were cal- 
culated via bootstrapping. In addition, in Figure[5]we in- 
dicate the rms field-to-field variations between the three 
sp atially distinct CO MBO-17 fields. As also pointed out 
in I Wolf et al.l l)2003j) . the field-to-field variations domi- 
n ate the error in the COMBO-17 jj^"*'' determinations. 

IBell et aP lj2003bj) point out that uncertainties in the 
absolute calibration and relative calibration of the SDSS 
and Johnson zeropoints can lead to < 10% errors in the 
derived rest-frame magnitudes and colors of galaxies. To 
account for this, we add a 10% error in quadrature with 
the formal errors for both the COMBO-17 and SDSS 
luminosity densities. These are the errors presented in 
Table [Hand Figure El 

In Figure we also plot the jy^^ value of luminous 
L BGs determined by in tegrating the luminosity function 
of IShaolev et^ pOOH) to L'^''=''''. A direct comparison 
between our sample and theirs is not entirely straight- 
forward because the LBGs represent a specific class of 
non-obscured, star forming galaxies at high redshift, se- 
lected by their rest-frame far UV light. Nonetheless, our 
jrcst determination at z — 2.8 is slightly higher than their 
determination at z = 3, indicating either that the HDF-S 
is overdense with respect to the area surveyed by Shapley 
et al. or that we may have galaxies in our sample which 
are not present in the ground-based LBG sample. 

D03 have also measured the luminosity density in the 
rest- frame _B-band but, because they do not give their 
luminosity function parameters except for their lowest 
redshift bin, it is not possible to overplot their luminosity 
density integrated down to our IV limit. 

4.LL The Evolution 0/ 

We find progressively stronger luminosity evolution 
from the V to the [/-band: whereas the evolution is quite 
weak in U, it is very strong in U . The f^^*^ in our high- 
est redshift bin is a factor of 1.9 ± 0.4, 2.9 ± 0.6, and 
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Table 1. Rest-Frame Optical Luminosity Density and Integrated 
Color 

logjfj-t logig'^' logj{?=t iU-B)rest {B-V)rest 

[hvo L0,tfMpc-3] [h7o Lq,bMpc-3] [hvo Lq^vMpc-^] 
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Note. — j'^^^ and rest-frame colors calculated for galaxies with 1V> 1.4 X 10^" hr^.^ Lq^v 
^SDSS 

bCOMBO-17 
<=FIRES 



4.9 ± 1.0 higher than the z = 0.1 value m V, B, and U 
respectively. To address the effect of cosmic variance on 
the measured evolution in j'^"*^* r ely on the clusterin g 
analysis developed for our sample in lDaddi et alJ l)2003j) . 
Using the correlation length estimated at 2 < z < 4, 
To — 5.5 hiQQ Mpc, we calculated the expected 1-sigma 
fluctuations in the number density of objects in our two 
highest redshift bins. Because for our high-z samples the 
poissonian errors are almost identical to the bootstrap 
errors, we can use the errors in the number density as a 
good proxy to the errors in j]^^*". The inclusion of the 
effects of clustering would increase the bootstrap errors 
on the luminosity density by a factor of, at most, 1.75 
downwards and 2.8 upwards. This implies that the in- 
ferred evolution is still robust even in the face of the 
measured clustering. The COMBO-17 data appear to 
have a slightly steeper rise towards higher redshift than 
our data, however there are two effects to remember at 
this point. First, our lowest redshift point averages over 
all redshifts z < 1.6, in which case we are in reasonably 
good agreement with what one would predict from the 
average of the SDSS and COMBO-17 data. Second, our 
data may simply have an offset in density with respect 
to the local measurements. Such an offset affects the 
values of but as we will show in i)4.2l it does not 

strongly affect the global color estimates. Nonetheless, 
given the general increase with j"^^^ towards higher red- 
shifts, we fit the changing j^^^^ with a power law of the 
form j]^^^{z) — j]^{0) * (1 + z)^. These curves are over- 
plotted in Figure 13 and the best fit parameters in sets of 
(ir'(O), /3) are (5.96 x 10^ 1.41), (6.84 x 10^ 0.93), and 
(8.42 X 10^, 0.52) in the U, B, and V bands respectively, 
where j^x^^{0) has units of hyp LgMpc^^. At the same 
time, it is important to remember that our power law 
fit is likely an oversimplification of the true evolution in 

■rest 
■^^ ■ . 

The increase in jJv with decreasing cosmic time can 
be modeled as a simple brightening of L, . Performing a 
test similar to that performed in ROl, we determine the 
increase in L^y with respect to L|f needed to match 
the observed increase in j'^^^^^ from z = 0.1to2.41<z< 
3.2, assuming the SDSS Schechter function parameters. 



To convert between the Schechter function para meters 
in the SDSS bands and those in the iBesseU l)1990D filters 
we transformed the L^'JJ.^^ values to the Bessell V filter 

using the {V — "'^ r) color, where the color was derived 
from the total luminosity densities in the indicated bands 
(as given in B03). We then applied the appropriate AB 
to Vega correction tabulated in Bessell (1990). Because 
the difference in Ag// is small between the two filters in 
each of these colors, the shifts between the systems are 
less than 5%. The luminosity density in the y-band at 
2.41 < z < 3.2 is j^""* = 1.53±0.26 x 10** hyo Lq Mpc-^. 
Using the 1^-band Schechter function parameters for our 
cosmology, (t>f^^^ = 5.11 x 10~^ hfg Mpc"^, a^°^^ = 
-1.05, and L^DSS ^ 2.53 x 10^° hf^f Lq, we can match 

the increase in jy^^ if L,.y brightens by a factor of 1.7 
out to 2.41 < z < 3.2. 

4.2. The Cosmic Color 

Using our measures of j"''' we estimated the cosmic 
rest-frame color of all the visible stars which lie in galax- 
ies with 1V> 1.4 X 10^*^ hyQ^ L0. We derived the mean 
cosmic {U — B)rest and {B — V)rest by using the j™*** 
estimates from the previous section with the appropri- 
ate magnitude zeropoints. The measured colors for the 
FIRES data, the COMBO-17 data, and the SDSS data 
are given in TableQ] Emission line corrected {U — B)rest 
colors may be calculated by applying Equation ^ to the 
values in the table. For the FIRES and COMBO-17 data, 
uncertainty estimates are derived from the same boot- 
strapping simulation used in S I4.1I In this case, however, 
the COMBO-17 and SDSS errorbars do not include an 
extra component from errors in the transformation to 
rest-frame luminosities, since these transformation errors 
may be correlated in a non-trivial way. 

The bluing with increasing redshift which could have 
been inferred from FigureElis seen explicitly in FigureOl 
The color change towards higher redshift occurs more 
smoothly than the evolution in f^^^, with our FIRES 
data meshing nicely with the COMBO-17 data. It is im- 
mediately apparent that the rms field-to-field errors for 
the COMBO-17 data arc much less than the observed 
trend in color, in contrast to Figure El This explicitly 
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Fig. 2. — The rest-frame optical luminosity density vs. cosmic 
age and redshift from galaxies with Kl°^^ < 25 and lV>Ly"'°'''\ 
For comparison we plot j'^*^*^* determinations from other surveys 
down to our L^'''^' limits. The squares are those from our data, 
the triangles are from the Combo-17 survey (Wolf ct al. 2003), the 
circle is that at z: = 0.1 from the SDSS (BOS), and the pentagon 
is that from ISTiaplev et al.l 120011) . The dotted crrorbars on the 
COMBO-17 data indicate the rms field-to-field variation derived 
from the three spatially distinct COMBO-17 fields. The solid line 
is a power law fit to the FIRES, COMBO-17, and SDSS data of 
the form i^'==*(z) = f^^'HO) * (1 + z)'^ ■ 
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shows that the integrated color is much less sensitive than 
jrcst .j-p field-to-field variations, even when such variations 
may dominate the error in the luminosity density. The 
COMBO-17 data at z < 0.4 are also redder than the 
local SDSS data, possibly owing to the very small cen- 
tral apertures used to measure colors in the COMBO-17 
survey. The colors in the COMBO-17 data were mea- 
sured with the package MPIAPHOT using using the peak 
surface brightness in images smoothed all to identical 
seeing (I'.'S). Such small apertures were chosen to mea- 
sure very precise colors, not to obtain global color es- 
timates. Because of color gradients, these small aper- 



tures can overestimate the global colors in nearby well 
resolved galaxies, while providing more accurate global 
color estimates for the more distant objects. Following 
the estimates of this bias provided bv .Bell et al. (2^0Q3y , 
we increased the color errorbars on the blue side to 0.1 
for the z < 0.4 COMBO-17 data. It is encouraging to 
see that the color evolution is roughly consistent with a 
rather simple model and that it is much smoother than 
the luminosity density evolution, which is more strongly 
affected by cosmic variance. 

We interpreted the color evolution as being primarily 
driven by a decrease in the stellar age with increasing 
redshift. Applying the {B — V)rest dependent e mission 
line corrections inferred from local samples (See i)3.2.1|l . 
we see that the effect of the emission lines on the color 
is much less than the magnitude of the observed trend. 
We can also interpret this change in color as a change in 
mean cosmic M/L* with redshift. In this picture, which 
is true for a variety of monotonic SFHs and extinctions, 
the points at high redshift have lower M/L* than those 
at low redshift. At the same time, however, the evolution 
in jy with redshift is quite weak. Taken together this 
would imply that the stellar mass density is also de- 
creasing with increasing redshift. We will quantify this 
in g3| 

To show how our mean cosmic (U — B)rest and 
{B — V)rest colors Compare to those of morphologically 
normal nearby galaxies, we overplot them in Figure |J 
on the locus of nearby galaxies from Larson & Tinsle^ 
ljl978j) . The integrated colors, at all redshifts, l i e very 
close to the local track, which iLarson fc TinslevI itTOT^) 
demonstrated is easily reproduceable with simple mono- 
tonically declining SFHs and which is preserved in the 
presence of modest amounts of reddening, which moves 
galaxies roughly parallel to the locus. In fact, correct- 
ing our data for emission lines moved them even closer 
to the local track. While we have suggested that M/L* 
decreases with decreasing color, if we wish to actually 
quantify the M/L* evolution from our data we must 
first attempt to find a set of models which can match 
our observed colors and which we will later use to con- 
vert between the color and M/Ly. We overplot in Fig- 
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ure ^ two model tracks corresponding to an exponen- 
tial SFH with T = 6 Gyr and with E(B - V) = 0, 
0.15, and 0.35 (assuming a iCalzetti et al.1 1)20001) red- 
dening law). These tracks were calcul a ted u sing the 
2000 version of the iBriizim l fc Chariot (1993) models 
and have Z = Zq and a ISalpete r (,.1955) IMF with a 
mass range of O.1-12OM0. Other exponentially declin- 
ing models and even a constant star forming track all 
yield similar colors to the r = 6 Gyr track. The mea- 
sured cosmic colors at z < 1.6 are fairly well approx- 
imated either of the reddening models. At z > 1.6, 
however, only the E{B — V) = 0.35 track can repro- 
duce the data. This high extinction is in contrast to 
the results of POl andlShaolcv ct al. (2001) who found a 
mean reddeni ng for LBGs of E(B - V) ~ 0.15 . SY98 
and Irhompso n. Wevmann. fc Storrie-LombardU |2001), 
however, measured extinctions on this order for galax- 
ies in the HDF-N, although the mean ext inction from 
iThomnson. Wevmann. fc Storrie-Lombardil pOOl) was 
lower at z > 2. The amount of reddening in our sample 
is one of the largest uncertainty in deriving the M/L* 
values, nonetheless, our choice of a high extinction is the 
only allowable possibility given the integrated colors of 
our high redshift data. 

Although this figure demonstrates that the measured 
colors can be matched, at some age, by this simple 
E{B — V) ^ 0.35 model, we must nevertheless inves- 
tigate whether the evolution of our model colors are also 
compatible with the evolution in the measured colors. 
This is shown by the track in Figure |21 We have tried 
different combinations of r, E{B — V), and Zstart, but 
have not been able to find a model which fits the data 
well at all r edshif ts. The parameterized SFR(z) curve of 
iCole et all 1)200 1|) also provided a poor fit to the data. 
Given the large range of possible parameters, our data 
may not be sufficient to well constrain the SFH. 

4.3. Estimating M/Ly and The Stellar Mass Density 

In this subsection we describe the use of our mean cos- 
mic color estimates to derive the mean cosmic M/Ly and 
the evolution in p». 

The main strength of considering the luminosity den- 
sity and integrated colors of the galaxy population, as 
opposed to those of individual galaxies, lies in the simple 
and robust ways in which these global values can be mod- 
eled. When attempting to derive the SFHs and stellar 
masses of individual high-redshift galaxies, the state-of- 
the-art models for the broad-band colors only consider 
stellar populations with at most two separate compo- 
nents (SY98; POl; Shapley et al. 2001). Using their stel- 
lar population synthesis modeling, Shaplcv ct al. (2001) 
proposes a model in which LBGs likely have smooth 
SFHs. On the other hand, SY98 concluded that 
they may only be seeing the most recent episode of 
star formation and that LBGs may indeed have burst- 
ing SFHs. This same idea was suppor ted by POl 
and lFerguson. D ickinso n, fc Papovich! f2002') using much 
deeper NIR data. When using similar simple SFHs to 
model the cosmic average of the galaxy population, a 
more self-consistent approach is possible. While individ- 
ual galaxies may, and probably do, have complex SFHs, 
the mean SFH of all galaxies is much smoother than that 
of individual ones. 

Encouraged by the general agreement between the 
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Fig. 3. — The evolution of the cosmic color plotted against red- 
shift and cosmic time for our data in addition to data from other 
^ < 1 surveys. The squares are those from our d ata, the triangles 
are from the Combo-17 survey iWoIf et alJI2no3) . and the circle is 
that at z = 0.1 from the SDSS (BOS). The open symbols indicate 
the empirical emission line correction to the integrated colors de- 
rived using the spectroscopic and photometric data from the NFGS 
iJansen ct al. 2000b). The dotted errorbars on the COMBO-17 
data indicate the field-to-field variation. Note that the integrated 
rest-frame color is much more stable than J^"''' against field-to-field 
variations. The COMBO-17 data point at 2; = 0.3 has been given a 
color errorbar of 0.1 in the blueward direction and an open symbol 
to reflect the possible systematic biases resulting from their very 
small central apertures. We also overplot a model with an expo- 
nentially declining SF H with r = %Gyr, Ej B - V) = 0.35, and 
Z3tart = 4.0 assuming a lCalzetti et al.l 120001) extinction law. 



measured colors and the simple models, we attempted to 
use this model to constrain the stellar mass-to-light ra- 
tio M/Ly in the rest-frame F-band, by taking advantage 
of the relation between color and logioM/L* found by 
iBell fc de Joii^ ()2001)) . For monotonic SFHs, the scatter 
of this relation remains small in the presence of modest 
variations in the reddening and metallicity because these 
effects run roughly parallel to the mean relation. Using 
the T — 6 Gyr exponentially declining model, we plot in 
Figure the relation between {U — V)rest and M/Ly 
for the E{B - F) = 0, 0.15, and 0.35 models. It is seen 
that dust extinction moves objects roughly parallel to the 
model tracks, reddening their colors, but making them 
dimmer as well and hence increasing M/Ly. Nonethe- 
less, extinction uncertainties are a major contributor to 
our errors in the determination of M/Ly. We chose to 
derive M/Ly from the {U — V)rest color instead of from 
the {B — V)rest color because at blue colors, where 
our high redshift points he, {B — V)rest derived M/Ly 
values are much more sensitive to the exact value of the 
extinction. This behavior likely stems from the fact that 
the {U — V)rest color spans the Balmer/4000A break 
and hence is more age sensitive than (B — V)rest- At 
the same time, while {U — B)rest colors are even less 
sensitive to extinction than [U — V)resti they are more 
susceptible to the effects of bursts. 

We constructed our relation using a ISalpeteil I^T955^ 
IMF^^. The adoption of a different IMF would simply 

We do not attempt to model an evolving IMF although evi- 
dence for a top-heavy IMF at high redshifts has been presented by 
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Fig. 4.— The {U - B)reat vs. [B - V)rest at z =1.12, 
2.01, and 2.8 of all the relatively unobscured stars in galaxies with 
1V> 1.4 X IQlO h^Q^ L0. The thick solid black Une is the local 
relation derived by Larson & Tinslcv 1 1978) from nearby morpho- 
logically normal galaxies. The large symbols are identical to those 
in Figure 121 For clarity we do not plot the field-to- field errorbars 
for the COMBO-17 data. The s mall solid points are the colors of 
nearby galaxies from the NFGS tjansen et a l. 2000a), which have 
been corrected for emission lines. The small crosses are the NFGS 
galaxies which harbor AGN. The thin tracks correspond to an ex- 
ponentially declinin g SFH with a ti mescale of 6 Gyr. The tracks 
were created using a'Salpeter" J1953) IMF and the 2000 version of 
the Bruzual & Chariot 1 1993) models. The dotted track has no ex- 
tinction, the dashed track has been reddened by E{B — V) = 0.15, 
and the th in solid track has be en reddened by E(B — V) = 0.35, 
using the ICalzetti et alj 120001) extinction law. The black arrow 
indicates the reddening vector applied to the solid model track. 
The emission line corrected data lie very close the track defined 
by observations of local galaxies and the agreement with the mod- 
els demonstrates that simple SFHs can be used to reproduce the 
integrated colors from massive galaxies at all redshifts. 



change the zeropoint of this curve, leaving the relative 
M/L* as a function of color, however, unchanged. As 
discussed in i)4.2l this model does not fit the redshift 
evolution of the cosmic color very well. Nonetheless, the 
impact on our M/Ly estimates should not be very large, 
since most smooth SFHs occupy very similar positions in 
the M/h*y vs. U-V plane. 

This relation breaks down in the presence of more com- 
plex SFHs. We demonstrate this in Figure where we 
plot the T = 6 Gyr track and a second track whose SFH 
is comprised of a 50 Myr burst at t = 0, a gap of 2 Gyr, 
and a constant SFR rate for 1 Gyr thereafter, where 50% 
of the mass is formed in the burst. It is obvious from this 
figure that using a smooth model will cause errors in the 
M/Ly estimate if the galaxy has a SFR which has an 
early peak in the SFH. At blue colors, such a early burst 
of SFR will cause an underestimate of M/Ly, a result 
similar to that of POl and DOS. At red colors, however, 
M/Ly would be overestimated with the exact systematic 
offset as a function of color depending strongly on the de- 
tailed SFH, i.e. the fraction of mass formed in the burst, 
the length of the gap, and the final age of the stellar 
population. 

The models show that the method may over- or under- 
estimate the true stellar mass-to- light ratio if the galaxies 
have complex SFHs. It is important to quantify the er- 

IFerguson. Dickinson, fc Papovichl 120021) 




U-V 



Fig. 5. — The relation between {U — V) and M/L^ for a model 
track with an exponential timescale of 6 Gyr. The dotted line is 
for a model with E{B — V) = 0, the dashed line for a model with 
E{B — V) = 0.15, and the solid line is for a model reddened by 
E{B — V) = 0.35 (using a Calzetti extinction law), which we adopt 
for our M/Ly conversions. The vertical solid arrows indicate the 
colors of the three FIRES data points, the vertical dotted arrow 
indicates the color of the SDSS data, and the diagonal solid arrow 
indicates the vector used to redden the E(B — V) = model to 
E{B — V) = 0.35. The labels above the vertical arrows correspond 
to the redshifts of the FIRES and SDSS data. 



rors on the global M/Ly based on the mean {U — V)rest 
color and how these errors compare to those when deter- 
mining the global M/Ly value from individual M/Ly 
estimates. To make this comparison we constructed a 
model whose SFH consist of a set of 10 Myr duration 
bursts separated by 90 Myr gaps. We drew galaxies at 
random from this model by randomly varying the phase 
and age of the burst sequence, where the maximum age 
was 4 Gyr. Next we estimated the total mass-to-light 
ratios of this sample by two different methods; first we 
determined the M/Ly for the galaxies individually as- 
suming the simple relation between color and mass-to- 
light ratio, and we took the luminosity weighted mean of 
the individual estimates to obtain the total M/Ly. This 
point is indicated by a large square in Figure Hand over- 
estimates the total M/Ly by ~ 70%. Next we first add 
the light of all the galaxies in both U and V , then use the 
simple relation between color and M/Ly to convert the 
integrated {U — V)rest into a mass-to-light ratio. This 
method overestimates M/Ly by much less, ~ 35%. This 
comparison shows clearly that it is best to estimate the 
mass using the integrated light. This is not unexpected; 
the star formation history of the universe as a whole is 
more regular than the star formation history of individ- 
ual galaxies. If enough galaxies are averaged, the mean 
star formation history is naturally fairly smooth. 

Using the relationship between color and M/Ly we 
convert our {U — V)rest and jy measurements to stel- 
lar mass density estimates . The resulting values are 
plotted vs. cosmic time in Figure |H1 We have included 
points for the SDSS survey created in an analogous way 
to those from this work, i.e. using the M/Ly derived 
from the rest-frame color and multiplying it by jy^^^ for 
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Fig. 6. — The effect of an early burst of star formation on 
the relation between {U — V) and M/L^. The relation between 
(U — V) and M/hy for a model track with an exponential timescale 
of 6 Gyr is show by the solid line. We also show a track for a SFH 
which includes a 50 Myr burst at 4 = followed by a gap of 2 
Gyr and then a constant SFR rate for 1 Gyr thereafter, where the 
fraction of mass formed in the burst is 0.5. The track continues for 
a total time of 4.5 Gyr. The dots are placed at 100 Myr intervals 
and the dotted section of the line indicates the very rapid transition 
in color caused by the onset of the second period of star formation. 
Both tracks have the same extinction. 
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Fig. 8. — The build-up of the stellar mass density as a func- 
tion of redshift. The solid points are for galaxies with 1V> 1.4 X 
10^'' h^Q Lq and were derived by applying the E(B — V) = 0.35 
relation in Figure |5] to the {U — V)rest colors and j?'^''^* mea- 
surements from the FIRES (solid squares) and SDSS data (solid 
circle). The y-axis scale on the left side corresponds to the p, 
values for these points. The hollow points show the total stellar 
mass density measurements from the one-component models in the 
HDF-N (D03; hollow stars; calculated assuming solar metallicity), 
the CFRS (Brinchmann & Ellis 2000; hollow circles), and the 2dF- 
GRS -I- 2MASS (Cole et al. 2001; hollow hexagon). The dotted 
errorbars on the DOS points reflect the systematic mass uncertain- 
ties resulting from metallicity and SFH changes. The y-axis scale 
on the right hand side corresponds to the p* estimates for these 
points. The relative scaling of the two axes was adjusted so that 
our SDSS p, estimate was at the same height as the total p» es- 
tima te of Cole e t al. T he solid curve is an integral of the SFR(2) 
from lTlole et alJ 1200 Ifl which has been fit to extinction corrected 
data at 21 < 4. 
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Fig. 7. — A comparison of different measures of the global 
M/Ly for a mock catalog of galaxies with bursting SFHs. The 
solid line represents the relation between {U — V) and M/L^ for 
a model track with an exponential timescale of 6 Gyr. The black 
dots show the true M/L's of the model starbursting galaxies, as 
described in the text; the open circle shows the true luminosity 
weighted Mtot/Ltot of the mock galaxies. The square shows the 
luminosity weighted M/L^ derived by applying the simple model 
to the individual galaxies - in this case, the mean Mtot/Ltot is 
overestimated by 70%. The triangle is the Mtot/Ltot derived from 
the luminosity weighted mean color (or [U — V)tot) of the model 
galaxies. It overestimates Mtot/Ltot by only 35%. 



all galaxies with lV>Ly'''^^'\ The estimates are listed 
in Table [5] We have derived the statistical errorbars on 
the yO* estimates by creating a Monte-Carlo simulation 
where we allowed jy^' and {U — V)rest (and hence 
M/Ly) to vary within their errorbars. We then took the 
resulting distribution of values and determined the 
68% confidence limits. As an estimate of our systematic 
uncertainties corresponding to the method we also deter- 
mined M/Ly from the {U — B)rest and (B — V)rest 
data using an identical relation as for the {U — V)rest 
to conversion. The {U - B).rest derived M/L^, 

values were different from the {U — V)rest values by 
a factor of 1.02, 0.80, 0.95, and 1.12 for the z =0.1, 
1.12, 2.01, and 2.8 redshift bins respectively. Likewise 
the {B — V)rest determined M/Ly values changed by 
a factor of 0.99, 1.11, 1.15, and 0.80 with respect to the 
{U — V)rest values. While the {U — V)rest values are 
very similar to those derived from the other colors, the 
{U — V)rest color is less susceptible to dust uncertainties 
than the {B — V)rest data and less susceptible to the 
effects of bursts than the {U — B)rest data. 

The derived mass density rises monotonically by a fac- 
tor of ~ 10 all the way to z 0.1, with our low redshift 
point meshing nicely with the local SDSS point. 

5. DISCUSSION 
5.1. Comparison with other Work 
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Table 2. M/L:j^ and Stellar, Mass Density Estimates 



2 log M/L^ log p* 



0.1 ±0.1=' 0.541^-^3 8.491; 

1.12t?-48b 0.131^6 8-14t: 

2.0lt»i;'' -0.421^8 7.48t! 

2.80lE;-^0b _o.701«-i^ 7.49t! 
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Figure |S1 shows a consistent picture of the build-up of 
stellar mass, both for the luminous galaxies and the to- 
tal galaxy population. It is remarkable that the results 
from different authors appear to agree well given that the 
methods to derive the densities were different and that 
the fields are very small. 

We compared our results to the total mass estimates 
of other authors in Figure |S1 In doing this we must re- 
member, because of our L^'^''^'^'^ cut, that we are missing 
significant amounts of light, and hence, mass. Assuming 
the SDSS luminosity function parameters, we lose 46% 
of the light at 2 = 0. Atz = 2.8, however, we inferred a 
brightening of y by a factor of 1.7, implying that we 
go further down the luminosity function at high redshift, 
sampling a larger fraction of the total starlight. If we ap- 
ply this brightening to the SDSS L,,y we miss 30% of the 
light below our luminosity threshold at z = 2.8. Hence, 
the fraction of the total starlight contained in our sample 
is rather stable as a function of redshift. To graphically 
compare our data to other authors we have scaled the 
two different axes in Figure |S1 so that our derivation of 
the SPS S Pi, is at t he sam e height of the total esti- 
mate of iCole et al.l 1)200 1|) . At z < 1 we compa r ed ou r 
mass estimates to those of iBrinchmann fc EllisI l)200(lD . 
Following DOS, we have corrected their published points 
to total masses by correcting them upwards by 20% to 
account for their mass incompleteness. The fraction of 
the total stars formed at z < 1 agrees well between our 
data and that of SDSS and Brinchmann & Ellis. At 
z > 0.5, we compared our results to those of D03. DOS 
calculated the total mass density, using the integrated lu- 
minosity density in the rest-frame B-band coupled with 
M/L measurements of individual galaxies. The fractions 
of the total stars formed in our sample (60%, 13%, and 
9% at z =1.12, 2.01, and 2.8) are almost twice as high as 
those of DOS. The results, however, are consistent within 
the errors. 

We explored whether field-to-field variations may play 
a role in the discrepancy between the two datasets. DOS 
studied the HDF-N, which has far fewer "red" galax- 
ies than HDF-S (e.g., Labbe et al. 2003, Franx et al, 
2003). If we omit the J — K selected galaxies found by 
iFranx et all l|2003D in the HDF-S, the formal M/L^, de- 
creases to 45% and 43% of the total values and the mass 
density decreases to 57% and 56% of the total values in 
the z = 2.01 and z = 2.8 bins respectively, bringing our 
data into better agreement with that from DOS. This re- 
inforces the earlier suggestion bv lFranx et al. ( 2003) that 
the J — K selected galaxies contribute significantly to the 



stellar mass budget. 

The errors in both determinations are dominated by 
systematic uncertainties, although our method should be 
less sensitive to bursts than that of DOS as it uses the 
light integrated integr ated over the galaxy population. 

We note that Fontan aet al.l l)200 3^ have also measured 
the stellar mass density in the HDF-S using a catalog 
derived from data in common with our own. They find a 
similar, although slightly smaller evolution in the stellar 
mass density, consistent with our result to within the 
uncertainties. 

5.2. Comparison with SFR(z) 

We can compare the derived stellar mass to the mass 
expected from determinations of the SFR as a function of 
redshift. We use the curve by Cole et al. (2001), who fit- 
ted the observed SFR as determined from various sources 
at z < 4. To obtain the curve in Figure |H1 we integrated 
the SFR(z) curve taking into account the time depen- 
dent stellar mass loss derived from the 2000 version of the 
J3ruzual fc Chariot (199S) population synthesis models. 

We calculated a reduced of 4.S when comparing 
all the data to the model . If, however, we omit the 
IBrinchmann fc Elli^ l)2000(l data, the reduced de- 
creases to 1.8, although the results at z > 2 lie sys- 
tematically below the curve. This suggests that some 
systematic errors may play a role, or that the curve is 
not quite correct. The following errors can influence our 
mass density determinations: 

-Dusty, evolved populations: it is assumed that the 
dust is mixed in a simple way with the stars, leading 
to a Calzetti extinction curve. If the dust is distributed 
differently, e.g., by having a very extincted underlying 
evolved population, or by having a larger R value, the 
current assumptions lead to a systematic underestimate 
of the mass. If an underlying, extincted evolved pop- 
ulation exists, it would naturally explain the fact that 
the ages of the Lyman-break galaxies are much younger 
than expected (e.g., POl, Ferguson et al. 2002). There 
may also be galaxies which contribute significantly to the 
mass density but are so heavily extincted that they are 
undetected, even with our very deep iiTs-band data. If 
such objects are also actively forming stars, they may be 
detectable with deep submillimeter observations or with 
rest-frame NIR observations from SIRTF. 

-Cosmic variance: the two fields which have been stud- 
ied are very small. Alth ough we use a consistent estimate 
of clustering from Daddi et alJ l)2003D . red galaxies make 
up a large fraction of the mass density in our highest red- 
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shift bins. Since red galaxies have a very high measured 
clustering from z ~ 1 (e.g., Daddi et al. 2000, Mc Carthy 
et al. 2001) up to possibly z - 3 l|Daddi et alJl2003|) . 
large uncertainties remain. 

-Evolving Initial Mass Function: the light which we see 
is mostly coming from the most massive stars present, 
whereas the stellar mass is dominated by low mass stars. 
Changes in the IMF would immediately lead to different 
mass estimates but if the IMF everywhere is identical (as 
we assume), then the relative masses should be robust. If 
the IMF evolves with redshift, however, systematic errors 
in the mass estimate will occur. 

-A steep galaxy mass function at high redshift: if much 
of the UV light which is used to measure the SFR at high 
redshifts comes from small galaxies which would fall be- 
low our rest-frame luminosity threshold then we may be 
missing significant amounts of stellar mass. Even the 
mass estimates of D03, which were obtained by integrat- 
ing the luminosity function, are very sensitive to the faint 
end extrapolation in their highest redshift bin. 

5.3. The Build-up of the Stellar Mass 

The primary goal of measuring the stellar mass den- 
sity is to determine how rapidly the universe assembled 
its stars. At z ~ 2 — 3, our results indicate that the 
universe only contained ^ 10% of the current stellar 
mass, regardless of whether we refer only to galaxies at 
1V> 1.4 X 10^" h^Q^ L0 or whether we use the total mass 
estimates of other authors. The galaxy population in 
the HDF-S was rich and diverse at z > 2, but even so it 
was far from finished in its build-up of stellar mass. By 
z 1, however, the total mass density had increased to 
roughly half its local value, indicating that the epoch of 
1 < z < 2 was an important period in the stellar mass 
build-up of the universe. 

A successful model of galaxy formation must not only 
explain our global results, but also reconcile them with 
the observed properties of individual galaxies at all red- 
shifts. For example, a population of galaxies at z ^ 
1 — 1.5 has been discovered (the so called extremely red 
objects or EROs), roughly half of which can be fit with 
formation redshifts higher than 2.4 ijCimatti et alJl2b02(l 
and nearly passive stellar evolution thereafter. Our re- 
sults, which show that the universe contained only ^ 10% 
as many stars at z 2 — 3 as today would seem to indi- 
cate that any population of galaxies which formed most 
of its mass at z > 2 can at most contribute ~ 10% of the 
present day stellar mass density. At z ~ 1 — 1.5, where 
the EROs reside, the universe had assembled roughly half 
of its current stars. Therefore, this would imply that the 
old EROs contribute about ^ 20% of the mass budget 
at their epoch. Likewise, it should be true that a large 
fraction of the stellar mass at low redshift should reside 
in objects with mass weighted stellar ages correspond- 
ing to a formation redshift of 1 < z < 2. In support of 
this,^ogE ct al. (2002) recently have shown that 40% 
of the local luminosity density at 0.7/iTO, and perhaps 
~ 50% of the stellar mass comes from centrally con- 
centrated, high surface brightness g alaxies whic h have 
red col ors. In ag reement with the 'Ho gg et alJ ( 2002 ) 
resufts, iBell et all (|2003ar and Ikauffm ann et alJ lj2003fi 
also found that ~ 50 — 75% of the local stella r mass 
density resides in early type galaxies. iHogg et aP 1)20021) 
suggest that their red galaxies would have been formed 



at z > 1, fully consistent with our results for the rapid 
mass growth of the universe during this period. 

6. SUMMARY & CONCLUSIONS 

In this paper we presented the globally averaged rest- 
frame optical properties of a Kg-hand selected sample 
of galaxies with z < 3.2 in the HDF-S. Using our very 
deep 0.3 — 2.2^m, seven band photometry taken as part 
of the FIRE Survey we estimated accurate photomet- 
ric redshifts and rest-frame luminosities for all galaxies 
with i^*°AB < 25 and used these luminosity estimates to 
measure the rest- frame optical luminosity density j^^^^, 
the globally averaged rest-frame optical color, and the 
stellar mass density for all galaxies at z < 3.2 with 
1V> 1.4 X 10^° hyp Lq. By selecting galaxies in the 
rest-frame y-band, we selected them in a way much less 
biased by star formation and dust than the traditional 
selection in the rest-frame UV and much closer to a se- 
lection by stellar mass. 

We have shown that f^^^^ in all three bands rises out 
to z - 3 by factors of 4.9 ± 1.0, 2.9 ± 0.6, and 1.9 ± 0.4 in 
the U, B, and y-bands respectively. Modeling this in- 
crease in j^^^^ as an increase in L» of the local luminosity 
function, we derive that must have brightened by a 
factor of 1.7 in the rest-frame T^-band. 

Using our j^^^^ estimates we calculate the {U — B)rest 
and [B — V)rest colors of all the visible stars in galaxies 
with 1V> 1.4 X 10^° hf(f Lq. Using the COMBO-17 data 
we have shown that the mean color is much less sensitive 
to density fluctuations and field-to-field variations than 
either f^^*^ or p*. Because of their stability, integrated 
color measurements are ideal for constraining galaxy evo- 
lution models. The luminosity weighted mean colors lie 
close to the locus of morphologically i iorma l local galaxy 
colors defined by iLarson fc TinslevI (IT971 . The mean 
colors monotonically bluen with increasing redshift by 
0.55 and 0.46 magnitudes in (U — B)rest and {B — V)rest 
respectively out to z 3. We interpret this color change 
primarily as a change in the mean stellar age. The joint 
colors can be roughly matched by simple SFH models if 
modest amounts of reddening {E{B — V)< 0.35) are ap- 
plied. In detail, the redshift dependence of {U — B)rest 
and (B — V)rest cannot be matched exactly by the sim- 
ple models, assuming a constant reddening and constant 
metallicity. However, we show that the models can still 
be used, even in the face of these small disagreements, to 
robustly predict the stellar mass-to-light ratios M/Ly of 
the integrated cosmic stellar population implied by our 
mean rest-frame colors. Variations in the metallicity does 
not strongly affect this relation and it holds for a variety 
of smooth SFHs. Even the IMF only affects the nor- 
malization of this relation, not its slope, assuming that 
the IMF everywhere is the same. The reddening, which 
moves objects roughly along this relation is, however, a 
large source of uncertainty. Using these M/Ly estimates 
coupled with our measurements, we derive the stel- 
lar mass density . These globally averaged estimates of 
the mass density are more reliable than those obtained 
from the mean of individual galaxies determined using 
smooth SFHs, primarily because the cosmic mean SFH 
is plausibly much better approximated as being smooth, 
whereas the SFHs of individual galaxies are almost defi- 
nitely not. 

The stellar mass density, p*, increases monotonically 
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with increasing cosmic time to come into good agreement 
with the other measured values at z < 1 with a factor of 
~ 10 increase from z ~ 3 to the present day. Within the 
ra ndom uncertainties, our results agree w ell with those 
of Dickinson. Papovich. Ferguso n, fc Biidavar i (2003) in 
the HDF-N although our estimates are systematically 
higher than in the HDF-N. Taken together, the HDF-N 
and HDF-S paint a picture in which only ~ 5 — 15% of 
the present day stellar mass was formed by z 2. By 
2 ~ 1, however, the stellar mass density had increased to 
~ 50% of its present value, implying that a large fraction 
of the stellar mass in the universe today was assembled 
at 1 < z < 2. Our estimates slightly underpredict 
the stellar mass density derived from the integral of the 
SFR(z) curve at z > 2. A resolution of the small ap- 
parent discrepancy between different fields, and between 
the predictions from optical observations will in part re- 
quire deeper NIR data, to probe further down the mass 
function, and wider fields with multiple pointings to con- 
trol the effects of cosmic variance. In addition, large 
amounts of follow-up optical/NIR spectroscopy are re- 
quired to help control systematic effects in the Zphot esti- 
mates. The 25 square arcminute MS1054-03 data taken 
as part of FIRES and the ACS/ISAAC GOODS obser- 
vations of the CDF-S region will be very helpful for such 
studies. Observations with SIRTF will also improve the 
situation by accessing the rest-frame NIR, where obscu- 
ration by dust becomes much less important. Finally, 
systematics in the M/L* estimates may exist because of 
a lack of constraint on the faint end slope of the stellar 
IMF. 

We still have to reconcile global measurements of the 
galaxy population with what we know about the ages 
and SFHs of individual galaxies. Our globally deter- 
mined quantities are quite stable and may serve as robust 
constraints on theoretical models, which must correctly 
model the global build-up of stellar mass in addition to 
matching the detailed properties of the galaxy popula- 
tion. 



Planck-Institute for Astrophysics (MPA), New Mexico 
State University, University of Pittsburgh, Princeton 
University, the United States Naval Observatory, and the 
University of Washington. 



GR would like to thank Jarle Brinchmann and Frank 
van den Bosch for useful discussions in the process of 
writing this paper. Christian Wolf for providing addi- 
tional COMBO-17 data products, and Eric Bell and 
Marcin Sawicki for giving comments on an earlier ver- 
sion of the paper. GR would also like to acknowledge 
the generous travel support of the Lorentz center and 
the Leids Kerkhoven-Bosscha Fonds, and the financial 
support of Sonderforschungsbereich 375. 

Funding for the creation and distribution of the SDSS 
Archive has been provided by the Alfred P. Sloan Foun- 
dation, the Participating Institutions, the National Aero- 
nautics and Space Administration, the National Sci- 
ence Foundation, the U.S. Department of Energy, the 
Japanese Monbukagakusho, and the Max Planck Soci- 
ety. The SDSS Web site is http://www.sdss.org/ 

The SDSS is managed by the Astrophysical Re- 
search Consortium (ARC) for the Participating Institu- 
tions. The Participating Institutions are The Univer- 
sity of Chicago, Fermilab, the Institute for Advanced 
Study, the Japan Participation Group, The Johns Hop- 
kins University, Los Alamos National Laboratory, the 
Max-Planck-Institute for Astronomy (MPIA), the Max- 



Rudnick et al. 



15 



APPENDIX 

DERIVATION OF Zp^ot UNCERTAINTY 

Given a set of formal flux errors, one way to broaden the redshift confidence interval without degrading the accuracy 
(as noticed in ROl) is to lower the absolute of every (z) curve without changing its shape (or the location of the 
minimum). By scaling up all the flux errors by a constant factor, we can retain the relative weights of the points in the 
without changing the best fit redshift and SED, but we do enlarge the redshift interval over which the templates 
can satisfactorily fit the flux points. Since we believe the disagreement between Zspec and Zphot is due to our finite and 
incomplete template set, this factor should reflect the degree of template mismatch in our sample, i.e., the degree by 
which our models fail to fit the flux points. To estimate this factor we first compute the fractional difference between 
the model and the data A^^ for the j**^ galaxy in the i'^ filter, 

A.. = ~J''' ^ (Al) 

where are the predicted fluxes of the best-fit template combination and J'^^* are our actual data. For each galaxy 
we calculated 



A, = 



\ 



-, JVfilt 

V A2 . (A2) 

A^mt-1^ 



where we have ignored the contributions of the [/-band. While the [/-band is important in finding breaks in the SEDs, 
the exact shapes of the templates are poorly constrained blueward of the rest-frame [/-band and the t/-data often 
deviates significantly from the best-fit model fluxes. 

To determine the mean deviation of all of the flux points from the model Adcv we then averaged over all galaxies 
in our complete FIRES sample with -ft'*°AB < 22 (for which the systematic Zphot errors should dominate over those 
resulting from photometric errors) to obtain 

A,e. = ^El^l- (A3) 

We find Adcv — 0.08, which includes both random and systematic deviations from the model. We modified the 
Monte-Carlo simulation of ROl by calculating, for each object j. 



l^i=2 \Sfi , . , 

^ ^fiit-1 ^ ' 



again excluding the [/-band. We then scaled the flux errors, for each object, using the following criteria: 

: 
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> > 
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/s\ r • (A5) 

^/^Adcv 

The photometric redshift error probability distribution is computed using the ^//'s. Note that this procedure will not 
modify the Zphot errors of the objects with low S/N where the Zphot errors are dominated by the formal photometric 
errors. The resulting probability distribution is highly non-Gaussian and using it we calculate the upper and lower 
68% confidence limits on the redshift ^^^qj and z'^f^^ respectively. As a single number which encodes the total range 

of acceptable Zphot's, we define Szphot = 0.5 * {Zpl^^ - Zp'^^t). 

Figure 6 from LOS shows the comparison of Zphot to Zgpec- For these bright galaxies, it is remarkable that our new 
photometric redshift errorbars come so close to predicting the difference between ZpUot and the true value. Some 
galaxies have large Szphot values even when the local x^ minimum is well defined because there is another minimum 
of comparable depth that is contained in the 68% redshift confidence limits. There are galaxies with Szphot < 0.05. 
Some of these are bright low redshift galaxies with large rest-frame optical breaks, which presumably place a strong 
constraint on the allowed redshift. Many of these galaxies, however, are faint and the Szphot is unrealistically low. 
Even though these faint galaxies have < they still can have high S/N in the -6450 or Vqqq bandpasses and 

hence have steep x^ curves and small inferred redshift uncertainties. In addition, many of these galaxies have Zphot 
> 2 and very blue continuum longward of Lya. The imposed sharp discontinuity in the template SEDs at the onset 
of HI absorption causes a very narrow minimum in the x^ (^) curve, and hence a small Szj,fiot, but likely differs from 
the true shape of the discontinuity because we use the mean opacity values ofjMadau (1995), neglecting its variance 
among different lines of sight. 

It is difficult to develop a scheme for measuring realistic photometric redshift uncertainties over all regimes. The 
Szphot estimate derives the Zphot uncertainties individually for each object, but can underpredict the uncertainties in 
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some cases. Compared to the technique of ROl however, a method based completely on the Monte-Carlo technique 
is preferable because it has a straightforwardly computed redshift probability function. This trait is desirable for 
estimating the errors in the rest-frame luminosities and colors and for this reason we will use Szphot as our uncertainty 
estimate in this paper. 

REST-FRAME PHOTOMETRIC SYSTEM 

To d efine the rest-frame U, B, and V fluxes we use the filter transmission curves and zeropoints tabulated in lBessell 
(|199(?), specifically his UX, B, and V filters. The Bessell zeropoints are given as magnitude offsets with respect to a 
source which has constant fi, and AB = 0. The AB magnitude is defined as 

AB, = -2.5 * logio(/,) - 48.58 (Bl) 

where (f^) is the flux f^iv) observed through a filter T^v) and in units of ergs s~^cm^^Hz^^. Given the zeropoint 
offset ZPi, for a given filter, the Vega magnitude rrii, is then 

m, = AB, - ZP, = -2.5 * logio(/.) - 48.58 - ZP,. (B2) 

All of our obs erved fluxes and rest-frame template fluxes are expressed in f\ . To obtain rest-frame magnitudes in the 
iBesselll ()1990j) system, we must calculate the conversion from f\ to fi, for the redshifted rest-frame filter set. The flux 
density of an SED with /a (A) integrated through a given fllter with transmission curve T(A) is 

^^'^ - jT'iX')dX' ^^^> 
^T'{u')dv' 

h{X')T'{X!)dX' = I Uv')T'{v')dv' (B5) 



or 



Since 



we can convert to (/y) through 

("HA>.{^ (B6, 

and use (/j^) to calculate the apparent rest-frame Vega magnitude through the redshifted filter via Eg. IB2I 

ESTIMATING REST-FRAME LUMINOSITIES 

We derive for any given redshift, the relation between the apparent AB magnitude of a galaxy through a 
redshifted rest-frame filter, its observed fluxes {f\i,obs) in the different filters i, and the colors of the spectral templates. 
At redshift z, the rest-frame filter with effective wavelength Xrest has been shifted to an observed wavelength 

A. = A^e.t X (1 + z) (CI) 

and we define the adjacent observed bandpasses with effective wavelengths A; and A;i which satisfy 

A; < A, < A;,.. (C2) 

We now define 

Cobs = mobs,\i — 'mobs,Xh (C3) 

where mobs,Xi and rnobs,Xh are the AB magnitudes which correspond to the fluxes (/a,.o6s) and (/a^.o^s) respectively. 
We then shift each template in wavelength to the redshift z and compute, 

where mtempi,Xi and mtempi,Xh are the AB magnitudes through the A; and Xt observed bandpasses (including the 
atmospheric and instrument throughputs). We sort the templates by their Ctempi values, Ctempi,a, Ctempi.b, etc., and 
flnd the two templates such that 

Ctempl,a Cobs < Ctempl,b- (^5) 

We then define for the a'^ template 

Cx„ (C6) 

where mtempi.x^ is the apparent AB magnitude of the redshifted a'^ template through the redshifted Xrest filter. 
We point out that because our computations always involve colors, they are not dependent on the actual template 
normalization (which cancels out in the difference). Taking our observed color Cobs and the templates with adjacent 
"observed" colors Ctempi,a and Ctempi,b, we can interpolate between CA,,z,a and C\i^z,b 

mobs,Xi — = Cxi^zM + (Cobs — Ctempl.a) X ( „ ■^^'^'^ Xi,z,a \ ^^^^ 

\^tenipl^b ^templ,a / 
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and solve for . 

When Cobs lies outside the range of the Ctempi^s, we simply take the two nearest templates in observed Ctempi space 
and extrapolate Eg. IC7l to compute m\^. 

Equation IC7I has the feature that « mobs,\i when — A/ (and hence when C\^^z.a and C\^^z.b ~ 0). While 
this method still assumes that the templates are reasonably good approximations to the true shape of the SEDs it has 
the advantage that it does not rely on exact agreement. Galaxies whose observed colors fall outside the range of the 
templates can also be easily flagged. A final advantage of this method is that the uncertainty in toa^ can be readily 
calculated from the errors in the observed fluxes. 

From mx^ , we compute the rest-frame luminosity by applying the K-correction and converting to luminosity units 

^ io-0.4(mA, - M0,A,.., - ^Pa_J ^ ( J^) ' x (1 + z)-' x h^' (C8) 
L© \10pcJ 

where ^^^^^ is the absolute magnitude of the sun in the Xrest filter (Mq u = +5.66, Mq b = +5.47, and 
M0^v = +4.82 in Vega magnitudes; Cox 2000), ZP\^^^^ is the zeropoint in that filter (as in Eq. 1621 but expressed at 
A and not at z/), and Dl is the distance modulus in parsecs. Following ROl, we correct this luminosity by the ratio 
of the _fir*°* flux to the modified isophotal aperture flux (see LOS). This adjustment factor, which accounts both for 
the larger size of the total aperture and the aperture correction, changes with apparent magnitude and it ranges from 
1.23, at 20 <i^*°AB < 24, to 1.69, at 24 <Kl°lj^ < 25, and it has an RMS dispersions of 0.17 and 0.49 in the two 
magnitude bins respectively. 

The uncertainty in the derived L'^''''* has contributions both from the observational flux errors and from the redshift 
uncertainty, which causes A^ to move with respect to the observed filters. The first effect is estimated by propagating 
the observed flux errors through Ea. lC7l As an example, object 531 at Zp^ot =2.20 has if^°AB — 24.91 and signal-to- 
noise in the if^-band of 8.99 and 5.43 in our modified isophotal and total apertures respectively. The resultant error 
in IV purely from fiux errors is then 26%. At i^'°AB ~ ^4, the typical signal-to-noise in the i^T^-band increases to w 13 
and w 6.3 in our modified isophotal and total apertures respectively and the error 17^^^ decreases accordingly. 

To account for the redshift dependent error in the calculated luminosity, we use the Monte-Carlo simulation described 
first in ROl and updated in [0 For each Monte-Carlo iteration we calculate the rest-frame luminosities and determine 
the 68% confidence limits of the resulting distribution. The 68% confidence limits in L^f^' can be highly asymmetric, 
just as for Zphot- For objects with if^°AB ^ ^5 we find that the contributions to the total L^'^''* error budget are 
dominated by the redshift errors rather than by the flux errors. 
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